In vitro and in silico evaluations of actinomycin X2and actinomycin D as potent anti-tuberculosis agents

Background Multidrug-resistant tuberculosis (MDR-TB) is one of the world’s most devastating contagious diseases and is caused by the MDR-Mycobacterium tuberculosis (MDR-Mtb) bacteria. It is therefore essential to identify novel anti-TB drug candidates and target proteins to treat MDR-TB. Here, in vitro and in silico studies were used to investigate the anti-TB potential of two newly sourced actinomycins, actinomycin-X2 (act-X2) and actinomycin-D (act-D), from the Streptomyces smyrnaeus strain UKAQ_23 (isolated from the Jubail industrial city of Saudi Arabia). Methods The anti-TB activity of the isolated actinomycins was assessed in vitro using the Mtb H37Ra, Mycobacterium bovis (BCG), and Mtb H37Rv bacterial strains, using the Microplate Alamar Blue Assay (MABA) method. In silico molecular docking studies were conducted using sixteen anti-TB drug target proteins using the AutoDock Vina 1.1.2 tool. The molecular dynamics (MD) simulations for both actinomycins were then performed with the most suitable target proteins, using the GROningen MAchine For Chemical Simulations (GROMACS) simulation software (GROMACS 2020.4), with the Chemistry at HARvard Macromolecular Mechanics 36m (CHARMM36m) forcefield for proteins and the CHARMM General Force Field (CGenFF) for ligands. Results In vitro results for the Mtb H37Ra, BCG, and Mtb H37Rv strains showed that act-X2 had minimum inhibitory concentration (MIC) values of 1.56 ± 0.0, 1.56 ± 0.0, and 2.64 ± 0.07 µg/mL and act-D had MIC values of 1.56 ± 0.0, 1.56 ± 0.0, and 1.80 ± 0.24 µg/mL respectively. The in silico molecular docking results showed that protein kinase PknB was the preferred target for both actinomycins, while KasA and pantothenate synthetase were the least preferred targets for act-X2and act-D respectively. The molecular dynamics (MD) results demonstrated that act-X2 and act-D remained stable inside the binding region of PknB throughout the simulation period. The MM/GBSA (Molecular Mechanics/Generalized Born Surface Area) binding energy calculations showed that act-X2 was more potent than act-D. Conclusion In conclusion, our results suggest that both actinomycins X2 and D are highly potent anti-TB drug candidates. We show that act-X2is better able to antagonistically interact with the protein kinase PknB target than act-D, and thus has more potential as a new anti-TB drug candidate.


INTRODUCTION
an essential tool in medical, developmental and molecular genetics as it is able to crosslink into duplex DNA, inhibit DNA-dependent RNA polymerase, and inhibit protein synthesis (Chen et al., 2012). In one study, act-D suppressed the coxsackievirus B3 and was proposed to treat AIDS (acquired immunodeficiency syndrome) due to its ability to block the minus-strand transfer mechanism in HIV-1 (Human Immunodeficiency Virus-1) (Chen et al., 2012).
We previously isolated, identified, and characterized a novel actinomycete strain, known as S. smyrnaeus strain UKAQ_23 (Fig. S1), which was obtained from a mangrove sediment sample in the Jubail industrial city of Saudi Arabia. We also produced, purified, identified, and characterized actinomycins X 2 and D, which are synthesized by the S. smyrneus strain UKAQ_23 (Figs. S2-S3) (Qureshi et al., 2021a). This current study is an extension of our previous work.
Here, we investigate the anti-TB potential of actinomycins X 2 and D (produced by S. smyrnaeus strain UKAQ_23) through in vitro and in silico studies, using 16 known anti-TB drug targets. Also, we show the mode of action for the anti-TB activity of act-X 2 and act-D against Mtb. This is the first study using molecular docking and dynamics studies to demonstrate the anti-TB potential of act-X 2 and act-D.

Test organisms
The Mtb H37Ra, BCG, and Mtb H37Rv were used as test organisms.

Isolation and identification of S. smyrnaeus strain UKAQ_23
The novel S. smyrnaeus strain UKAQ_23 was isolated from a mangrove sediment sample collected from the Jubail industrial city of Saudi Arabia (Fig. S1). The strain UKAQ_23 was identified as a novel S. smyrnaeus strain UKAQ_23 by 16S rRNA gene sequencing. The 16S rRNA gene sequences of UKAQ_23 were submitted to GenBank (NCBI) with accession number MG657032.1. A culture copy of S. smyrnaeus strain UKAQ_23 has also been deposited at the National Centre for Microbial Resource (NCMR), Pune, MS, India, under the accession number MCC 0192. Our previous publication provides detailed information on this novel strain UKAQ_23 (Qureshi et al., 2021a).

Production, extraction, purification, and identification of actinomycins X 2 and D
The novel strain UKAQ_23 was used to produce actinomycins using a modified ISP-4 fermentation medium and the solid-state fermentation method (Fig. S2). The actinomycins were extracted, purified, and identified as act-X 2 and act-D. Our previous publication provides more detailed information on act-X 2 and act-D (Qureshi et al., 2021b;Qureshi et al., 2021a).

In vitro Anti-TB activity of act-X 2 and act-D
The in vitro anti-TB activity of act-X 2 and act-D was conducted using the Microplate Alamar Blue Assay (MABA) method (Collins & Franzblau, 1997;Franzblau et al., 1998;Amin et al., 2008;Mendoza-Aguilar et al., 2012;Cho, Lee & Franzblau, 2015). The minimum inhibitory concentration (MIC) values for act-X 2 and act-D against Mtb H37Ra and BCG were determined using a single read-out assay (fluorescence), while the MIC values for act-X 2 and act-D against Mtb H37Rv were calculated using a dual read-out assay (optical density and fluorescence). The mean of the two MIC values (from both single and dual read-out assays) was used. Each test was performed in triplicate. The results are expressed using the mean ± standard deviations (SD).

Microplate Alamar Blue Assay (MABA) cell stock preparation
All solutions and media were sterilized and stored at 4 • C away from direct sunlight until use. The following materials were used in the preparation of MABA cell stocks: sterile 20% (v/v) tween 80, 50% (v/v) glycerol, Middlebrook 7H11 agar medium, Middlebrook 7H9 broth medium, phosphate buffer saline (PBS), PBS-Tween 80 (PBST), tryptic soy agar with 5% sheep blood. All procedures were carried out in a sterile environment. The inocula for the test organism's seed stocks were prepared by suspending colonies from 7H11 agar in 200 mL of 7H9 broth in a 500 mL Nephelo flask, and allowing them to grow until the log phase (40-60 Klett units) was attained; after which they were removed from the flask. The suspensions were then transferred to 50 mL conical centrifuge tubes and centrifuged at 4,000 × g for 10 min at 4 • C. After centrifugation, the supernatants were discarded, and one mL of sterile 1 x PBS was dispensed into each of 50 mL conical centrifuge tubes, followed by resuspension of pellets by pipetting. Each 50 mL tube was then filled with 30 mL of PBS, and pellets were resuspended in PBS by vortexing for 5 s. All tubes were centrifuged at 4,000 × g for 10 min at 4 • C. After centrifugation, the supernatants were discarded, and each 50 mL conical centrifuge tube was filled with one mL of PBS, followed by resuspension of pellets by pipetting. Aliquots containing 200 µL of suspended cells were transferred into sterile screwcap micro-tubes with a capacity of 1.5 mL. The 200 µL sample was then used to validate sterility on the blood agar plate and determine colony-forming units (CFU) on 6-well plates, following the serial dilution of 10 −4 -10 −7 . All tubes were labeled and stored at −80 • C in a storage box. The Klett units, tube numbers, sterility, CFU test results, date, and location in the freezer (−80 • C) were all recorded on a log sheet. The seed stocks were preserved for up to one year.

MABA protocol
The following materials were used in the MABA protocol: 20% (v/v) Tween 80, 7H12 medium, quality control (QC) standard anti-TB drugs, i.e., rifampin (0.02 µg/mL), isoniazid (0.46 µg/mL), linezolid (0.92 µg/mL), moxifloxacin (0.19 µg/mL), rifampicin (2 µg/mL), and streptomycin (6 µg/mL), sterile 96-well clear plates, a sonicator, blood agar plate, 6-well plates with five mL of 7H11 agar in each well, Alamar Blue ® , and a microplate fluorometer with 530 nm excitation and 590 nm emission filters. To prevent medium evaporation from the test wells during incubation, 200 µL of 7H12 medium was dispensed into the outer-perimeter wells of the 96-well plates. All wells in column 3 were filled with 200 µL of 7H12 medium, while all other wells were filled with 100 µL of 7H12 medium. The 2 µL of the listed QC drugs or test actinomycins, X 2 and D, were added in columns 1-3, rows B-F in a 96-well microtiter plate. Two-fold serial dilutions of act-X 2 and act-D were then prepared from columns 3 to 10, and at the end of dilution, 100 µL was discarded from column 10. The frozen cell stocks were thawed and sonicated briefly for 3-5 s at 20% power to achieve an even suspension. The cultures were diluted with the 7H12 medium until they reached a density between 5× 10 4 and 1× 10 5 CFU/mL. Then, 100 µL of culture was dispensed into wells from columns 2-11, rows B-G. The plates were incubated for seven days at 37 • C with 5% CO 2 and 95% humidity. The sterility of the cell stock was determined by streaking it onto a blood agar plate and incubating it at 37 • C for one week. Ten-fold serial dilutions were prepared up to 10 −7 , and 50 µL of each dilution was then inoculated into 6-well 7H11 agar plates and incubated for three weeks at 37 • C in a CO 2 incubator. If bacterial or fungal colonies were isolated on the blood agar or 7H11 plates within a week, the batch of cell stock was autoclaved and discarded. After three weeks, the CFU on the 7H11 agar plates was counted and reported. The 96-well plates were removed from the incubator on the seventh day of incubation. The Alamar Blue-Tween 80 solution was prepared by diluting Alamar Blue in an 8:5 ratio with 20% (v/v) Tween 80 (for one plate, 2 mL Alamar Blue and 1.25 mL 20% (v/v) Tween 80 solution were required). A 32.5 µL of Alamar Blue-Tween 80 was added to each well in the 96-well plates. Plates were incubated for 18-24 h at 37 • C, 5% CO 2 , and 95% humidity. The fluorescence was measured at 530 nm excitation and 590 nm emission. The MIC 90 was then determined using the following formula: The MIC 90 value is defined as the lowest concentration of a test compound that inhibits fluorescence by 90% compared to an untreated bacterial control.
Each protein was handled individually in the Biovia Discovery Studio Visualizer 2020 and MGLTools 1.5.6 programs, and standard procedures were followed for the preparation of receptors-including deletion of co-crystallized ligands, water molecules, and cofactors and addition of Gasteiger charges (Shushni, Azam & Lindequist, 2013;Hussain et al., 2020). Three-dimensional structures of the compounds act-X 2 and act-D were retrieved from the PubChem database. The ligand preparation module of MGL Tools 1.5.6 was used for merging all non-polar hydrogens, and to define rotatable bonds of the ligands.
AutoDock Vina 1.1.2 was used for molecular docking analysis using the default protocol (Trott & Olson, 2009). A grid box was built with a grid spacing of 1 Å at the center of respective co-crystallized ligands, with dimensions of 30 points in all directions for each receptor. Upon completing the docking, the best poses were selected from the top ten models from each target by examining their binding energy ( G binding , kcal/mol) and non-bonded interactions profile (Azam, Eid & Almutairi, 2021). Biovia Discovery Studio Visualizer 2020 and PyMol 1.7.4 were used for studying molecular interactions.

In silico molecular dynamics (MD) simulations
After the detailed docking analysis, the best complexes with lowest binding energies were selected. MD simulations of Protein Kinase PknB complexed with act-D and act-X 2 were then performed for 250 ns using the GROningen MAchine For Chemical Simulations (GROMACS) simulation software (GROMACS 2020.4) (Abraham et al., 2015), using the Chemistry at HARvard Macromolecular Mechanics 36m (CHARMM36m) (Huang et al., 2016) forcefield for proteins and the CHARMM General Force Field (CGenFF) for act-X 2 and act-D. In a truncated octahedral box, both protein-ligand systems were solvated with TIP3P water molecules. The minimum distance between protein systems and the edges of the simulation box was set to 10 Å, to efficiently meet the criteria for minimum image convention during MD simulation. The 0.15M KCl ionic concentration was introduced by adding 39 K + and 33 Cl − ions to the environment. The protonation states for His, Tyr, Lys, Arg, Asp, and Glu residues were evaluated at pH 7 and 0.15M of salinity using an H++ web server (version 3.2) and implemented after visual inspection (Anandakrishnan, Aguilar & Onufriev, 2012). The systems for act-D and act-X 2 contained a total of 39,707 and 39,787 atoms, respectively. The Chemistry at HARvard Macromolecular Mechanics Graphical User Interface (CHARMM-GUI) webserver was used to generate all input files (Jo et al., 2008;Lee et al., 2016).
Both systems were minimized for 5,000 steps using the Steepest Descent technique. Convergence was achieved under the force limit of 1,000 (kJ/mol/nm) to exclude any steric clashes. Later, both minimized systems were separately equilibrated at NVT (Canonical ensemble: where moles, N; volume, V; and temperature, T were conserved) and NPT (Isothermal-Isobaric ensemble: where moles, N; pressure, P; and temperature, T were conserved) ensembles for 100 ps (50,000 steps) and 1000 ps (1,000,000 steps), respectively, using time steps 0.2 and 0.1 fs, at 300 K to ensure a fully converged system for the production run (Qureshi et al., 2021a;Qureshi et al., 2021b).
The simulation runs for both systems were conducted at a constant temperature of 300 K and a pressure of 1 atm, or 1 bar (using an NPT ensemble), utilizing weak coupling velocity re-scaling (modified Berendsen thermostat) and Parrinello-Rahman algorithms, respectively. The relaxation times were set at τ T = 0.1 ps and τ P = 2.0 ps. Using the LINear Constraint Solver (LINCS) algorithm, all bond lengths involving hydrogen atoms were maintained stiffly at optimal bond lengths, with a time step of 2 fs. The non-bonded interactions were calculated using Van der Waals and Coulomb potentials. Interactions within a short-range cutoff of 12 Åwere calculated in each time step. The electrostatic interactions and forces in a homogeneous medium beyond the short-range cutoff were calculated using the Particle Mesh Ewald (PME) method. The Periodic Boundary Conditions (PBC) were applied in all x,y, and z directions. The production was run for 250 ns for each of the three complexes. The trajectory and energy data were recorded at every 10 ps (Qureshi et al., 2021b). The first 25 ns of the trajectories were excluded for comprehensive analysis, and the remaining 225 ns were utilized (except for MM/GBSA calculations, where the first 50 ns of the trajectories were excluded, and the remaining 200 ns were utilized). GRaphing, Advanced Computation and Exploration of data (Grace) was used to generate all plots (https://plasma-gate.weizmann.ac.il/Grace). The Molecular Mechanics/Generalized-Born Surface Area (MM/GBSA) protein-ligand binding energy was averaged over 200 frames calculated at the interval of 1ns over the last 200 ns of 250 ns production run (Miller et al., 2012). The following formula was used to calculate the binding free energy: Five additional MD simulations (replicates) for each of the two systems (act-X 2 and act-D) were also carried out to authenticate the ligand binding further. The starting conformation for each simulation was taken after every 50 ns for a 250 ns long trajectory (i.e., 50 ns, 100 ns, 150 ns, 200 ns and 250 ns). Before the final production run, each system was subjected to energy minimization and equilibration at both NVT and NPT ensembles. Each simulation was run for 50 ns, followed by MM/GBSA calculations. As previously mentioned, all other parameters and values remained the same.

In vitro anti-TB activity of act-X 2 and act-D
Act-X 2 and act-D exhibited substantial anti-TB activity against the tested mycobacterial strains. When tested against Mtb H37Ra, BCG, and Mtb H37Rv, act-X 2 exhibited MIC values of 1.56 ± 0.0, 1.56 ± 0.0, and 2.64 ± 0.07 µg/mL, while act-D showed MIC values of 1.56 ± 0.0, 1.56 ± 0.0, and 1.80 ± 0.24 µg/mL, respectively ( Fig. 1 and Table S2). Act-X 2 and act-D had similar MIC values against the Mtb H37Ra and BCG strains, although, the MIC value was lower when act-D was tested against the Mtb H37Rv strain. This indicates that actinomycins X 2 and D are more potent against the Mtb H37Ra and BCG strains than the Mtb H37Rv strain.

Statistical analysis
A statistically significant difference (p < 0.05) was observed between the anti-TB activity of act-X 2 and act-D, as determined by one-way ANOVA (analysis of variance); F (1, 2) = 22.471, p = 0.042 (Table S3).

In silico molecular docking studies
To gain insight into the potential interaction modes of act-X 2 and act-D, a molecular docking study was performed using several anti-TB drug targets (Timo et al., 2019;Baptista et al., 2021). The predicted binding energies of the act-X 2 and act-D compounds are presented in Table 1. Figures 2-3 and Figs. S4-S19 illustrate the binding energies of the docked molecules with the various anti-TB drug targets. The bar plot clearly shows variations in the affinities for the proteins studied. PknB plays a role in multiple mycobacterial biochemical pathways, and is crucial for in vitro growth and survival of the microbe in the host (Chawla et al., 2014). Of the 16 drug targets studied, PknB was the most promising target for both act-X 2 and act-D-with binding energies of −11.7 and −11.8 kcal/mol, respectively. In addition, both compounds, act-X 2 and act-D preferentially bound to DNA -indicating that they were equipotent in terms of docking, with a predicted binding affinity of −10.6 kcal/mol (Fig. S19).
KasA was the least favorable target protein for act-X 2 (binding energy: −5.7 kcal/mol), and pantothenate synthetase was the least favorable target for act-D (binding energy: −6.0 kcal/mol). Variable affinities were observed for all other targets ( Fig. S4 and Table 1).

In silico MD simulations
The trajectory analyses for act-X 2 and act-D in complex with protein kinase PknB were performed in terms of root mean square deviation (RMSD) of protein and ligand, root mean square fluctuations (RMSF), the radius of gyration (g (r) ), and the center-of-mass distance (CoM) between the protein and ligands (Fig. 4); their mean values are given in Table 2. The number of H-bonds, radial distribution function (RDF) of ligand around protein, the partial density of protein and ligands, and MM/GBSA (Miller et al., 2012) ligand-protein binding energies (Fig. 5) were also plotted. The values for various energy contributions are listed in Table 3. Figure 4A shows the RMSD plot for protein kinase (PknB) in complex with both actinomycins X 2 (red line) and D (black line). Both systems remained stable throughout the simulation time; although, the mean RMSD value for act-X 2 (2.85 ± 0.62 Å) was slightly higher than that of act-D (2.59 ± 0.31 Å). The RMSD value increased for act-X 2 between 175 to 200 ns, which was associated with the stronger movements of residues 157 to 180 forming a loop. Of note, this loop was absent from the crystal structure due to the higher beta-factor value. Overall, the RMSD values fluctuated within the range of ∼3 Å, indicating that both systems were stable and no major conformational changes were observed, other than the loop mentioned above. Figure 4B displays the RMSD plot for ligands actinomycins X 2 (red line) and D (black line). Both actinomycins X 2 and D remained highly stable within the binding pocket of the PknB protein, and exhibited mean RMSD values of 1.87 ± 0.27 Åand 1.77 ± 0.29 Å, respectively. Although act-X 2 displayed more conformational changes than act-D, the increase in stability was negligible-as reflected by its smaller standard deviation value.  Figure 4C shows the fluctuation in distance between the center-of-masses (CoM) of protein kinase (PknB) and actinomycins X 2 (red line) and D (black line) during the simulation time. Although act-X 2 had higher RMSD values, it remained completely intact with the PknB protein binding site, with a mean CoM distance of 15.96 ± 0.87 Å. Act-D partially exited the binding pocket during ∼80 ns to 175 ns of simulation time, with an overall mean CoM distance of 18.17 ± 2.79 Å, but still remained partially intact with the protein. This partial exit of act-D may explain why the ligand-induced local conformational changes in PknB were not as prominent as with act-X 2 . Figure 4D depicts the radius of gyration (g (r) ) plot for PknB in complex with actinomycins X 2 (red line) and D (black line) during the simulation time. Both protein systems remained highly compact, with mean values of 19.53 ± 0.27 Å and 19.68 ± 0.20 Å. They did not exhibit any significant global ligand-induced conformational changes. Figure 4E shows the root means square fluctuations for PknB in complex with actinomycins X 2 (red line) and D (black line). Similar patterns of fluctuation were observed for both protein systems. Both N-and C-terminals were flexible, as is commonly observed for globular proteins. Residues 13 to 23, which constitute a loop directly involved in ligand binding, showed ligand-induced fluctuations. These residues showed consistent interactions with act-D during its partial exit from the binding pocket of PknB, which resulted in slightly higher RMSF values compared to those observed for act-X 2 . Residues 157 to 180 form a loop when PknB is in complex with act-X 2, and this is comparatively more flexible than the complex with actinomycin D. This might be due to the reduction in contact time caused by the partial exit of act-D from the binding pocket of PknB, as mentioned earlier. Another significant fluctuation can be observed for residues 78 to 89, corresponding to another loop in the proximity of the N-terminal. Overall, the    Figures 5A-5B displays the number of hydrogen bonds formed during the simulation between PknB and ligands, actinomycins X 2 (red line) and D (black line). Act-X 2 formed consistent H-bonds with PknB throughout the simulation. In contrast, act-D only formed H-bonds during ∼80 ns to 175 ns of the simulation, when it partially exited the binding pocket of PknB. In fact, act-D seldom forms H-bonds with PknB. After 175 ns of simulation, the ligand formed a greater number of H-bonds with the protein, which may indicate that act-D and PknB have a slower binding rate than act-X 2 to reach the same efficacy; however, a detailed kinetic study would be needed to confirm this.

Table 3 MM/GBSA binding energy in kcal/mol for protein-ligand complexes (protein kinase PknB in complex with Act-D and Act-X2). (A)
MM/GBSA binding energy in kcal/mol for protein-ligand complexes (protein kinase PknB in complex with Act-D and Act-X2) averaged over 200 frames taken at the interval of 1ns over the last 200 ns of 250 ns production run. (B) MM/GBSA binding energy in kcal/mol for proteinligand complexes (protein kinase PknB in complex with Act-D and Act-X2) averaged over 50 frames taken at the interval of 1ns from 50 ns production run for 10 replicates. The different starting conformation for each simulation was taken from 250 ns long trajectory with an interval of 50 ns.  Figure 5C displays the radial distribution functions of ligands actinomycins X 2 (red line) and D (black line), around the PknB protein. The peak maxima for the radial distribution of act-X 2 lies close to 12.39 Å (higher intensities representing increased complex stability). A slight shoulder peak around 16.06 Å was observed, indicating a different distance for the protein-ligand complex. The peak maxima for the act-D-protein RDFs lie higher, close to 18.96 Å. However, an indistinguishable shoulder at a lower distance was also visible. Both RDFs show broad peaks, indicating that the ligands are flexible within the protein binding pocket. Figure 5D shows the partial density of PknB in complex with actinomycins X 2 (solid red line for protein and dashed line for ligand) and D (solid black line for protein and dashed line for ligand). No significant changes in the partial protein density were observed for either protein system. However, substantial changes in the distance between the peak maxima were noted for the partial density of ligands and their corresponding protein changes, indicating that the ligands may bind different regions of the binding site. Figure 5E shows the molecular mechanics/Generalized Born surface area (MM/GBSA) ligand-protein binding energies. The energy contributions are given in Table 3A. The binding performance was better between Act-X 2 and PknB ( G MM/GBSA = −26.10 ± 4.93 kcal/mol) compared with Act-D and PknB. Act-D initially showed a higher G, however, after 175 ns the binding energy was equivalent to that of act-X 2. The mean G MM/GBSA value was −19.04 ± 7.29 kcal/mol. This shows that, although act-D has the same potential to bind PknB as act-X 2 , its binding reaction rate is slower; however, further kinetics studies would be required to confirm this.
Van der Waals forces are essential for favorable ligand-protein binding, while polar solvation free energy is highly unfavorable. Coulombic energy and non-polar solvation free energy are also favorable for ligand-protein binding. To assess ligand binding, 10 different replicates of act-D and act-X 2 were tested to obtain 50 ns long trajectories for each simulation. These results complement our previous observations for both ligands (Table  3B).
Our in silico data supported the anti-TB activity we observed in vitro. The in silico data demonstrated that both actinomycins X 2 and D had potential anti-TB activity against PknB. Act-X 2 exhibited stronger binding affinity with faster kinetics, while act-D took longer to stabilize in the protein binding site with a slower reaction rate; although, further kinetics studies would be needed to confirm this.

DISCUSSION
MDR-TB is a highly contagious disease caused by MDR-Mtb. Due to the increasing prevalence of MDR in Mtb, treatment of MDR-TB is one of the most challenging issues affecting healthcare. Discovery of novel anti-TB drug candidates and targets is therefore essential to overcome these challenges. We therefore aimed to investigate the anti-TB efficacy of actinomycins, X 2 and D (obtained from the novel S. smyrnaeus strain UKAQ_23) in vitro, and to use in silico studies to further support these findings. The Mtb protein kinase B (PknB) has therapeutic potential-particularly as it plays a role in cell growth regulation and cell wall synthesis, and it induces morphological changes during bacterial cell division (Wehenkel et al., 2008;Khan, Kaur & Nandicoori, 2018). In addition, PknB exhibited the best binding affinity with both actinomycins X 2 and D. Therefore, PknB complexes with actinomycins X 2 and D were studied using molecular dynamics (MD) simulations to investigate their structure and dynamics, and to assess their potential for inhibiting PknB function. The simulation trajectories were further subjected to molecular mechanics/Generalized Born surface area (MM/GBSA) binding energy evaluations to find ligand-protein binding interactions and associated thermodynamics. Earlier, several studies exploited MD simulation for the development of potential antitubercular drug targets/molecules (Almeleebia et al., 2021;Burley et al., 2021). The in vitro anti-TB activity of act-X 2 and act-D was determined using the MABA method (the most widely adaptable method). We used three Mycobacterium strains: Mtb H37Ra, BCG, and Mtb H37Rv. Our results are consistent with previous studies of other strains that produce act-X 2 and act-D (i.e., Streptomyces sp. IMB094, Streptomyces sp. MS449, S.nasri YG62, S.elizabethii. II, S.padanus JAU4234, S.griseoruber, S.flavogriseus NJ-4, S.heliomycini, Streptomyces strain M7, S.hydrogenans IB310, Streptomyces MITKK-103, and Streptomyces sp. HUST012) (Kurosawa et al., 2006;Bird & Latif, 2007;El-Naggar, El-Assar & Abdul-Gawad, 2009;Praveen & Tripathi, 2009;Chen et al., 2012;Xiong et al., 2012;Khieu et al., 2015;Wei et al., 2017;Kulkarni et al., 2017;Wang et al., 2017a;Wang et al., 2017b;Sharma & Manhas, 2019). In addition, several previous reports demonstrated that act-X 2 and act-D had significant anti-TB activity (Praveen & Tripathi, 2009;Chen et al., 2012). Chen et al. (2012) reported the anti-TB activity of act-X 2 and act-D produced by Streptomyces sp. MS449 against BCG and Mtb H37Rv. They reported MIC values of 0.5 µg/mL for act-X 2 and 0.5 µg/mL for act-D against BCG; and 1.0 µg/mL for act-X 2 and 8.0 µg/mL for act-D against Mtb H37Rv (Chen et al., 2012). These findings are consistent with our results, demonstrating the substantial anti-TB activity of isolated act-X 2 and act-D against BCG and Mtb H37Rv. Another report demonstrated a MIC value of 0.78 µg/mL for the anti-TB activity of Py2 (act-D) against Mtb H37Rv, which is also consistent with our findings and further supports the anti-TB activity of act-X 2 and act-D against Mtb H37Rv (Praveen & Tripathi, 2009). Although these earlier studies demonstrated actinomycins X 2 and D showed anti-TB potential, they did not perform any in silico anti-TB investigations. Our study demonstrates that the act-X 2 and act-D, produced by the novel S. smyrnaeus strain UKAQ_23, exhibit substantial anti-TB activity against the Mycobacterium strains tested. In addition, we have demonstrated that isolated act-X 2 and act-D had comparable anti-TB efficiency against the Mtb H37Ra and BCG strains; although, act-D exhibited significantly higher anti-TB activity than act-X 2 against the Mtb H37Rv strain.

CONCLUSIONS
In conclusion, we suggest that both actinomycins X 2 and D, (newly sourced from the novel S. smyrnaeus strain UKAQ_23), are highly potent anti-TB drug candidates. Our results indicate that act-X 2 is more potent than act-D, and thus has more potential to be used as a new anti-TB drug candidate-by antagonistically interacting with the kinase PknB target protein.